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G.TANAHASHI 

Institute for Cosmic Ray Research, University of Tokyo 
Tanashi, Tokyo, 188 Japan 

ABSTRACT 

The distribution of the number of muons in EAS and the equi- 
intensity curves of EAS are analyzed on the basis of Monte Carlo 
simulation of various cosmic ray composition and the interaction 
models. Problems in the two best combined models are discussed. 

1 . Introduction 

Many analyses of EAS data have been reported to investigate the cosmic 
ray composition and their interaction. For example, the comprehensive 
work was reported by Gaisser et al(1) in 1978. After then, the equi- 
intensity data of Chacaltaya was revised (2) and the other data of Akeno 
experiment was added to it(3). Akeno experiment also gave a distribution 
of the number of muons (>1Gev) of EAS at fixed shower size with a good 
statistics (4) . In this report we make an analysis of these data with use 
of the Monte Carlo simulation described in 2. 

2. One dimensional simulation of EAS 

Among the parameters used in the interation model of the cosmic ray 
with the air nucleus, the collision m.f.p. is assumed to be 
/VA)/(1+0.07(logE^(Tev) + 1 ) 1 ,5 )gcm“ 2 where Ao=80 and 120 for proton and pion 
interaction respectively, and Eo is the cosmic ray energy. The value of 
A are shown in Fig.1 together with ^=290 .Eo (Tev)^ # ^^mb obtained in the 
Akeno experiment (5) . The assumptions of the other parameters are follow- 
ing. 

The leading particle carrys the energy Es fluctuated uniformly between 
0-Eo. The partition of the rest of the energy between the fragmentation 
and the central region is half to half. The energy spectrum of the pro- 
duced particles in the fragmentation is exp(-E/E*) where E* is 0.6Es and 
r Es for proton and pion incident respectively. In the central region the 
spectrum is assumed to be a plateau shaped one whose energy spread is 
shown in Fig. 2 for various models. Resultant multiplicity-energy relation 
, multiplicity distribution, scaling behavior and inclusive rapidity dis- 
* tribution are checked to fit the data of accelerator experiments. 

Produced particles are to be pions or kaons, and the proportion of the 
kaon production is to increase with the energy as 0 .OSlogE^Gev) . 

Interation models shown in Fig. 2 are used for the EAS calculation. 

They are SC(scaling), STD(standard, which is scaling + the extension of 
the central distribution in lO^ev to the higher energies), SCI /2 (seal ing 
+ increasing multiplicity with Ej / 2 in the central region above 10 1Zf ev), 
SCM (seal ing + particles nearly at rest in CM system above lO^ev) and CCM 
(no fragmentation particles and all particles nearly at rest in CM system 
above lO^ev). 

We have further two kinds, I and II, in each of the above models. I 
includes the generation of the high energy neutral pion as the leading 
particle, and the forward and backward symmetry of the produced particles 
in CM system. II includes no production of the energetic neutral pion as 
the leading , and the enhanced particle production in the backward three 
times more than the forward as the effect of target nucleus. The details 
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of these assumtions 
will be described els- 
where . 

The shower curves 
of electrons and muons 
of STD based on I, II 7 ~ 

are shown in Fig. 3. It 
is seen that the 
attenuation of particles 
between 1 000-1 600gcnT 2 6 - 

is nearly exponential 
and has almost the same w 
attenuation length (180- S 
190gcm“ 2 x and 190- w 5 - 
200gcm “ 2 in II) irres- | 
pective of the primary 55 
energies above 10 ^ev. S 
In the same figure, the 
shower curves of the 4 ~ 

constant cross section 
are also shown for the 
comparision. The reason 
that the attenuations — 

are almost the same Fig. 3 

comes mainly from the r — 

increasing cross sec- 
tions. The other result 
is that the muon con- 
tent in II is almost 3 7 _ 

times larger than I T s. 

In Fig. 4 the shower w 
curves of various h 

models are shown in & 
comparison. All of then 1 6 “ 
are based on model I. « 

It is found again the § 
attenuation byond 1000 
gcm “ 2 are almost para- 5 - 
llel. The difference 

among SC, STD and SCI/2 

are very small, and this Fig. 4 
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means that the energy distribution in the central region does not affect 
much the shower curve except the extreme case like SCM. The crosses in 
the figure show the iron shower curve (HSD) based on the superposition of 
STD. 

3. Nu distribution of fixed Ne 

The Np distributions of EAS 

of different primary component 
(mass number A=1,4,15 and 50) 
are calculated under the con- 
dition of fixed Ne assuming the 
primary integral energy spect- 
rum to be E“2. These Np dis- 
tributions include the error 00 

which comes into Nju in the p 

course of the data reduction g 
and also the error due to the | 
possible fluctuation of muon g 

lateral distribution. The dis- § 
tribution of A>4 are almost due 06 ioo 
to th&e errors. These N^u dis- 
tributions of each component 
are superposed to fit the ex- 
periment (5) as shown in Fig. 5. 

The conclusions obtained from 
this comparision are (i)low Nu 
tail requires the existence of 
proton showers more than Fig. 

several % (ii)the Nu distri- 
bution of the proton shower must be as broad as in STD, that is, the very 
high multiplicity models are avoided because of their narrow distribution 
of Np (iii)the avrage Np of each component must be spread over more than 
factor 3 for fixed Ne. 

4. Equi^ intensity curves of EAS 

Equi-intensity curves of a mixed composition are related to each shower 
curves as follows. Assume that the shower curves are energy independent 
in the required energy region, that is N(E,x)=E.N(x) where N(e,x) is the 
shower curve of energy E in depth of x. The intensity ratio of the heavy 
to the proton shower in fixed shower size is 
4(x)=(Ih/lp) N =(lh/Ip) E .(Nh(x)/Np(x)) Y 
where k=(lh/lp) E is the intensity ratio in fixed energy and V is the 
integral exponent of the primary energy spectrum. Putting Y =2, /KX) can 
be calculated with use^of the values of Nh(x), Np(x) as the funtion of k. 
Equij-intensity curves N(x) are obtained from 

N(x)=Np (x)/ (1+A(x))+Nh(x) . A(x)/(1+A(x)> 

A(x) and Tf(x) are calculated for the following combinations of the proton 
and the heavy showers. 


proton shower 

STD 

STD 

STD 

SCM 

heavy" shower 

HSD 

SCM 

CCM 

CCM 

notation 

HSD/STD 

SCM/ STD ■ 

CCM/STD 

CCM/SCM 


The N(x) are stiwn in Fig. 6 for k=0.1, 1.0 and 10 together with the ex- 
perimental data(2)(3). In general, the attenuation length between 800-1400 
gcm“2 are steeper than either prediction, but two cases (HSD/STD and k=10) 
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(CCM/SCM and k-0.1) barely fit the experiment as shown in the figure in 
dotted curves after the normalization of the intensity. SCM/STD and CCM/ 
STD do not fit because of too much difference of proton and heavy shower 
curves. This differenc makes the composite curves flat. 


5. Discussions and conclusion 

Both two selected cases satisfy ( the requirement (i) in 3. 
As for CCM/SCM, it is not 
clear whether (ii) (iii) are 



satisfied or not, but seems 
to be not. This model is 
introduced to approach the 
one which is proposed by 
Kakimoto et al(2) to explain 
the rapid development of EAS 
in high altitude observed by 
them. So the further pro- 
gress of this experiment 
will test the model. 

On the other hand, HSD/ 

STD satisfys all the require- 
ments (i) (ii) (iii) in 3. The 
problems of this model are 
the less content ofmuons in 
EAS and the slow development 
of the number of muons. The 
model needs about 2 times 
more muons . If we adopt 
model II instead of I, the 
number of muons increases 
but the attenuation length 
also increases apart fur- 
ther from the experiment. 

As for tfie .slow development 
of muon in this model, we 
can scarcely fit the experi- 
ment by taking the number 
of muons of iron showers in 
600 gcm“2 and of proton 

Showers in 1 100gcm~ 2 (Fig. 7) . If the N^i difference of the proton and the 
heavy nucleus shower is larger in some model, such a model is better in 
this problem. 
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